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We investigate the transfer of energy from large scales to small scales in fully developed forced 
three-dimensional MHD-turbulence by analyzing the results of direct numerical simulations in the 
absence of an externally imposed uniform magnetic field. Our results show that the transfer of 
kinetic energy from the large scales to kinetic energy at smaller scales, and the transfer of magnetic 
energy from the large scales to magnetic energy at smaller scales, are local, as is also found in the 
case of neutral fluids, and in a way that is compatible with Kolmogorov (1941) theory of turbulence. 
However, the transfer of energy from the velocity field to the magnetic field is a highly non-local 
process in Fourier space. Energy from the velocity field at large scales can be transfered directly 
into small scale magnetic fields without the participation of intermediate scales. Some implications 
of our results to MHD turbulence modeling are also discussed. 
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I. INTRODUCTION 

Most astrophysical and planetary systems, e.g. 
solar/stellar winds, accretion disks and interstellar 
medium, are in a turbulent state and coupled to mag- 
netic fields. Understanding and quantifying the statis- 
tical properties of magnctohydrodynamical (MHD) tur- 
bulence is crucial to explain many physicalprocesses in 
the cosmos, and in industrial flows as well p). Although 
the phenomenology of hydrodynamical (HD) turbulence 
is understood to some extent, and the theory has been 
able to make predictions like Kolmogorov's 4/5th law and 
the functional form of the energy spectrum in the inertial 
range, that have been well verified in experiments and nu- 
merical simulations, a similar statement cannot be made 
for MHD turbulence at the same level. In MHD flows, 
the two fields (velocity and magnetic) and two associ- 
ated energies involved in the dynamical processes allow 
for many possibilities for the energy to transfer between 
smaller or larger scales, making the dynamics more com- 
plex to address in both theory and modeling. 

We briefly describe some phenomenological aspects of 
HD turbulence to point out some of the difficulties usu- 
ally encountered when the formulation of HD turbulence 
is applied in the MHD case. To follow Kolmogorov (1941) 
theory Q (hereafter, K41), we need to assume a statis- 
tically isotropic and homogeneous flow in steady state in 
which the energy is cascading from eddies of scale I to 
smaller eddies, and so on until energy reaches the dis- 
sipation scales. Since we are considering a statistically 
steady state, the flux of energy to smaller scales has to 
be constant. We can further assume that the flux at some 
scale can depend only on the scale I and the amplitude 
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of the velocity field |u;| at this scale. This assumption 
is justified by the argument that larger eddies will only 
advect smaller eddies without significantly altering their 
scale and only when eddies of similar size interact do 
they produce a cascade. Therefore, only "local" inter- 
actions among the different scales control the cascade. 
Here we use the term "local" in terms of the different 
scales involved (i.e. scales of similar size) and not as 
locality in physical space. With these assumptions we 
obtain that the energy |u;| 2 at the scale I will cascade to 
smaller scales in a time i/|uj|, and since the energy cas- 
cade rate e is constant, we obtain e ~ | u.; | 3 /Z that implies 
|u;| ~ I 1 / 3 , which finally leads to the well verified K41 
spectrum dE/dk ~ k~ 5 ^ 3 to within small intermittency 
corrections. 

The assumptions of the HD theory of turbulence have 
been tested in the literature. Rcf. Q first tested the 
assumption of locality using direct numerical simulations 
(DNS) of 64 3 grid points. Their work has been followed 
by a number of authors with higher resolution simula- 
tions H H |1 Q, i, H E3- Refs. H [H GJ have 
also investigated the effect of long-range interactions and 
anisotropy induced by an anisotropic large scale flow. Al- 
though some issues still remain regarding the effect of 
long range interactions, the locality of the energy trans- 
fer has been confirmed. 

However, there are two important assumptions used in 
the HD case that are not necessarily true for the MHD 
case. First, the assumption of isotropy breaks down if 
an imposed uniform magnetic field is considered. We 
will not investigate such effects in the present work and 
will only consider flows with J b dx 3 = 0. The second 
assumption, that of locality of interactions among the 
different scales is what motivates our work. Unlike the 
HD case where the effect of larger eddies on smaller ones 
is the advection of the later ones (an effect that can be 
taken away by a Galilean transformation), in MHD the 
effect of a large scale fluctuation of the magnetic field 
cannot be so eliminated. Therefore, in MHD it is pos- 
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sible for small scales to interact directly with the large 
scales. If this is the case, we can not consider a "con- 
tiguous" transfer of energy in wave number space and 
cannot a priori follow the same arguments Kolmogorov 
used for HD turbulence. Therefore knowledge of the en- 
ergy transfer among different scales is important for the 
construction of any phenomenological model of turbu- 
lence. 

Present phenomenological models follow Kolmogorov 
like arguments that take into account the effect of the 
magnetic field. Iroshnikov ^3 and Kraichnan 0] pro- 
posed the first models to describe isotropic-MHD tur- 
bulence, predicting a spectrum of fc~ 3 / 2 (hereafter, IK). 
Goldreich and Shridar |l5l | proposed a new model for 
anisotropic MHD turbulence that takes into account the 
anisotropy introduced by a uniform magnetic field B , 
predicting a spectrum of k ± 5 ^ 3 , where k± refers to the 
direction perpendicular to Bo. Several models have 
been proposed that combine the two spectra (see e.g. 

El III Q), 

suggesting that the index of the energy 
spectrum is sensitive to the presence and intensity of 
Bo- Some aspects of non-locality of interactions are taken 
into account in the afore-mentioned models by consider- 
ing that large scale fluctuations of the magnetic field act 
as a uniform magnetic field to the smaller scales, and as 
a result they speed up or slow down the rate at which the 
energy is cascading. However, in these models, although 
non-local interactions are taken into account, the energy 
is transfered locally from one scale to a slightly smaller 
scale, like in Kolmogorov's HD turbulence model. 

The locality of the interactions and the energy transfer 
in MHD turbulence has been investigated through vari- 
ous closure models. The energy transfer has been stud- 
ied within the EDQNM closure model by E3 an d more 
recently by pol | where non-local interactions have been 
noted. Using field theoretical calculations the transfer of 
energy has been estimated by 0, 1^2, ■ As far as we 
know, the locality of the energy transfer in MHD has been 
investigated through three dimensional direct numerical 
simulations (DNS) only very recently [24| (see also j2f| 
for the two-dimensional case). These authors measured 
the transfer of energy between different scales and fields 
using free decaying MHD turbulence simulations with 
512 3 grid points. Their results showed that there is lo- 
cal transfer of energy between the same fields, while the 
transfers involving the two different fields showed a less 
local behavior, in the sense that a wider range of scales 
was involved in the interactions. 

In our work we use the results of DNS of mechanically 
forced MHD turbulence (unlike the free decaying case 
studied in |2^|) to study the locality of the energy trans- 
fer between different scales and fields. In all the cases 
studied we consider a mechanic external forcing that gen- 
erates a well defined large scale flow and small scale tur- 
bulent fluctuations. This is a regime of interest for several 
astrophysical and geophysical flows where magnetic fields 
are believed to be sustained against Ohmic dissipation by 
a dynamo process |26j . and the only external source of 



energy driving the system is mechanical (e.g. convection 
and rotation). There is an important difference between 
the case studied in [24| and the case considered in our 
work. In our case energy is forced through the velocity 
field and the system reaches a steady state with equipar- 
tition between the two fields. For this to happen there 
must be a non-zero flux for all times from the velocity 
field to the magnetic field. This is not necessarily true 
for the case of decaying turbulence and as our results 
show this significantly modifies the energy transfers from 
the velocity field to the magnetic field. 

In Sec. |H]we introduce the definitions of the transfer 
terms for MHD, and in Sec. IIIII we present the code we 
use for the numerical simulations as well as the results of 
the analysis. Finally, in Sec. II VI we summarize the main 
results of our work. 



II. THEORY AND DEFINITIONS 

The equations that describe the dynamics of an incom- 
pressible conducting fluid coupled to a magnetic field in 
the MHD approximation are given by: 

d t u + u • Vu = -Vp + b • Vb + vV 2 u + f, (1) 
d t h + u • Vb = b • Vu + ?yV 2 b, (2) 



V-u = 0, V-b = 0, (3) 

where u is the velocity field and b is the magnetic field. 
p is the (total) pressure and v and r\ are the viscosity 
and the magnetic diffusivity respectively. Here, f is the 
external force that drives the turbulence and the dynamo. 
The largest wavenumber of the Fourier transform of / is 
going to be denoted as kp and we are going to refer to 
kpl -1 as the forced scale. We are also going to define 
the viscous dissipation scale as fc" 1 = (e/^ 3 ) -1 / 4 and 
resistive dissipation scale as fc^ 1 = (e/r; 3 ) -1 / 4 where e is 
the energy dissipation rate. A large separation between 
the two scales (|k F | _1 ^> maxljfc" 1 , fc" 1 }) is required 
for the flow to reach a turbulent state. 

To investigate the transfer of energy among different 
scales of turbulence we use the Fourier transforms of the 
fields: 

u(x) =X>(kK kx , u(k) = ^JL y u(x)e- ikx dx 3 
and 

b(x) = £b(k) e ikx , b(k) = ^ J b(x)e- ikx dx 3 , 

where the domain is taken to be a triply periodic cube 
of size L = 2ir. We can now introduce the shell filter 
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decomposition: 

u ( x ) = u k(x), b(x) = bit (x) 



K 



K 



where 



%W = u(k)e 



K<\k\<K+l 

and similar for the field b 



'kx 



'kx 



K<\k\<K+l 

The fields and hx are therefore defined as the part 
of the velocity and magnetic field respectively, whose 
Fourier transform contains only wave numbers in the 
shell (K, K + 1] (hereafter called the shell K) and repre- 
sent "eddies" of scale K~ l . The evolution of the kinetic 
energy in a shell K, E U (K) = J u 2 K /2dx 3 is given by: 

d t E u {K) = /"^[-u K -(u-V)-u Q + u i f(b-V)-b t 



-v\S7u K \ + f • u K dx 3 



(4) 



and for the magnetic energy E b (K) = jh 2 K /2dx 3 we 
obtain 



E b(K) = |E[- b ^u-V)-b Q +b x .(b-V).u c 



-T?|Vb K | 2 dx 3 



(5) 



The above equations can be written in the more compact 
form: 

d t E u (K) = J2[Tuu(Q,K)+T bu (Q,K)]-vV u (K)+f(K), 
Q 

(6) 

d t E b {K) = YjFubiQ, K) + %b(Q, K)) - rfD h {K). (7) 
Q 

Here we have introduced the functions T UU (Q,K), 
T ub (Q, K), T bb (Q, K), and % U (Q, K) that express the en- 
ergy transfer between different fields and shells. 

T UU (Q, K) expresses the transfer rate of kinetic energy 
lying in the shell Q to kinetic energy lying in the shell K 
through the velocity advection term and is defined as: 



T UU {Q, K) = - J u K (u • V)u Q cZx 3 



(8) 



We similarly define 

T bb (Q,K) ee -j b K (u- V)b Q dx 3 , (9) 

which expresses the rate of energy transfer of magnetic 
energy lying in the shell Q to magnetic energy lying in 



the shell K through the magnetic advection term. The 
Lorentz force is responsible for the transfer of energy from 
the magnetic field to the velocity field. The resulting 
transfer rate is defined as: 



T bu (Q,K) ee J u K (b-V)b Q dx 3 



(10) 



Finally the term responsible for the stretching of the mag- 
netic field lines results in the transfer from kinetic energy 
to magnetic energy, given by: 



T ub (Q,K)= J b K (b.V)u Q dx 3 



(11) 



In summary, the functions T VW (Q, K) (for arbitrary fields 
v and w) represent the rate of transfer of energy from the 
field v (first index) in the shell Q (first argument), into 
energy of the field w (second index) in the shell K (second 
argument). If T VW (Q,K) > 0, then a positive amount of 
v— energy is transfered from the shell Q to w— energy in 
the shell K. If T VW (Q,K) < 0, then a negative amount 
of v— energy is transfered from the shell Q to w— energy 
in the shell K, or in other words, energy is transfered 
backwards from the shell K to the shell Q. 

In eqs. <|()I7[) we have also introduced two dissipation 
functions: the kinetic energy dissipation rate 



vV u (K) 



ivu^rrfx 3 , 



(12) 



and the magnetic energy dissipation rate 

TjVhiK) = v J \Vb K \ 2 dx 3 . (13) 



Finally, 



HK) 



f-u K d-x. 3 



(14) 



is the energy injection rate to the velocity field through 
the forcing term. 

Before presenting the results from numerical simula- 
tions, let us discuss some of the properties of the transfer 
functions. If T VW (Q, K) (where v,w can be either u or 
b) is expressing the rate of energy transfer from the field 
v in the shell Q to the field w in the shell K, then the 
following identity should hold 



T VW (Q,K) = -T WV (K,Q). 



(15) 



The interpretation of eq. H15|l is that the rate at which the 
shell Q is giving energy to the shell K must be equal to 
the rate the shell K is receiving energy from the shell Q. 
Eq. (|15|l can be easily shown to hold for all the transfer 
functions we defined (eqs. |8lllj ). It is this property that 
allows us to interpret the functions T UUl 7~ bu , T ub , and T bb 
as the energy transfer between different scales and fields. 

For a turbulent flow in a statistically steady state, 
equations © and J7J imply that: 

J2(Tuu(Q,K) + T bu (Q,K)) = (V U (K))-(F(K)), (16) 
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and 

^(T„ b (Q, K) + % b (Q, K)) = (V b (K)), (17) 
Q 

where (•} stands for a time average or an ensemble av- 
erage. For fixed K outside the forcing band, and in the 
limit of v, ?y — > 0, we have that 

Y J {Tuu{Q,K)+T hu {Q,K)) = Q (18) 
Q 

and 

Y,{Tub{Q,K)+T bb {Q,K)) = Q. (19) 
Q 

However, limited resolution will allow us to be in the 
regime where these last two equations hold only for a 
small range of wavenumbers. 

Finally we need to comment on the definitions of the 
various transfer functions we are using in this paper and 
the connection to the triad of wave numbers (k, p, q) that 
satisfy the relation k + p + q = (because of the convo- 
lution term resulting from the quadratic nonlinearities of 
the primitive equations); such triad is the basis for mode 
to mode interactions (see e.g. |27|). Our approach is 
equivalent to considering all triad interactions with the 
one wavenumbcr k G K and q £ Q and summing over 
all p satisfying k + p + q = in all shells, where p is the 
wave-number of the advecting field, and k and q are the 
wavenumbers of the modes energy is transfered to and 
from. Although the approach we are using gives us infor- 
mation on whether the energy is transfered locally or not, 
it cannot give definite conclusions on whether the inter- 
actions themselves are local. For example, even if energy 
is transfered locally from a wavenumbcr k to a wavenum- 
bcr q ~ k, the wavenumber p that is responsible for the 
transfer is not necessarily of the same order of magnitude 
as |k| and |q|. Ideally, one would investigate transfer 
terms of the form: T UU (K\P\Q) = J u^(up • V)ugdx 3 
that contain information about the third wave number 
involved in the interactions taking place. However the 
difficulty of manipulating data from high resolution runs 
and the difficulty of interpreting the results of transfer 
functions that depend on three arguments restricts us 
for the present time, to examine just the locality of the 
energy transfer. 

III. RESULTS 

To study the transfer of energy in MHD turbulence 
we use the turbulent steady state of several mechanically 
forced three dimensional MHD direct numerical simula- 
tions. The simulations and details of the code can be 
found in |28l . l29j| . The runs were performed in a triply pe- 
riodic domain with a resolution of 256 3 grid points, using 
a pseudo-spectral scheme with the 2/3-rule for dealiasing. 



TABLE I: Simulations. L is the integral length-scale of the 
flow, denned as L = J E u (k) dkj f E u (k)k~ 1 dk\ v is the kine- 
matic viscosity, and rj the magnetic difhisivity. The kinetic 
and magnetic Reynolds numbers Re and Km are based on L 
and the rms velocity, while the ratio of magnetic to kinetic 
energy Et/E u is the average in the turbulent steady state. 



Forcing 


L 


V 


V 


Re 


Rm 


Eb/E u 


ABC 


1.64 


2 x KT^ 


2 x irr a 


820 


820 


0.84 


TG 


1.35 


2 x 1(T 3 


5 x 1(T 3 


675 


270 


0.72 



The equations were evolved in time using a second order 
Runge-Kutta method. 

Turbulence was generated by two different types of 
forcing. In the first case a non-helical Taylor-Green 
force (hereafter referred as TG) was used itg(&o) = 
( sin(/coa;) cos(koy ) co s(fc n z), — cos(fc x) sin(fc y) cos(k z), 0) 
with k = 2 [231 . In the second case a helical 
ABC force was used fABc(ko) = {Bcos(k y) + 
Csin(fcoz), Ccos(fc z) + Asm(k x), Acos(kox) + 
Bsm(koy)) with fco = 2 poT| . All simulations were 
done with constant in time external force. First a 
hydrodynamic simulation was carried using each force, 
to reach a turbulent steady state. Both external forces 
generate a well defined large scale flow at \Kp\ ~ 3, 
and small scale turbulent fluctuations following to a 
good approximation a 5/3 Kolmogorov law. Then MHD 
simulations were carried, and a small magnetic field was 
amplified and sustained to equipartition by a dynamo 
process. The results in this paper are based on the 
saturated stage of the dynamo, which we will refer in 
the following as the MHD turbulent steady state. 

The transfers were calculated based on the definitions 
(|51- [TT|) . The transfer of energy during the early stages of 
the MHD simulations, when the magnetic energy is small 
and the velocity field is not modified by the Lorentz force 
(often referred to as the kinematic dynamo regime) are 
examined in a companion paper [30l | (hereafter referred 
as Paper II). Tablc^gives several relevant parameters for 
each run, and figure^Jshows the resulting energy spectra. 

Both simulations display a large scale magnetic field, 
although the spectrum of magnetic energy in the ABC 
simulation shows a stronger peak at k = 1. This peak 
is related with the dynamo a-effect and the inverse cas- 
cade of magnetic helicity. Details of this process will be 
discussed in Paper II. However, it is important to note 
that in the ABC simulation the large scale magnetic field 
is strongly helical, while in the TG simulation the mag- 
netic helicity is negligible. This large scale magnetic field 
is self-sustained by the turbulence. In both simulations, 
the net cross helicity (correlation between the velocity 
and the magnetic field) is small and can be neglected. 
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FIG. 1: Spectra of kinetic energy (solid line) and magnetic en- 
ergy (dashed line) of the ABC and Taylor Green runs, where 
the Taylor Green spectra have been shifted down by a factor 
of 20 for clarity. The Kolmogorov slope is showed as a ref- 
erence. Note that the magnetic Prandtl number Pm = vjr\ 
differs for the two runs. 
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FIG. 2: The transfer of energy T U u(Q,K) for the Taylor- 
Green run. The figure shows the rate that energy is transfered 
from the modes Q = 3, 10, 20, 30 to all the other modes K. 
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FIG. 3: The transfer of energy T UU (Q, K) for the ABC run. 
The figure shows the rate that energy is transfered from the 
modes Q — 3, 10, 20, 30 to all the other modes K. 



slightly smaller wavenumbers (negative T uu ) and give en- 
ergy to modes with slightly larger wave number (positive 
T uu ). The locality of the transfer is expressed from the 
fact that the transfer of energy from the modes in the 
shell Q to modes in shells K with K <C Q or K 3> Q 
is very small, and decreases fast with the separation of 
the two wave numbers. Finally, as the shell wavenumbcr 
K and Q is increased, there is a drop in the amplitude 
of the transfer. If the transfer functions were self-similar 
then an increase of the wave numbers K and Q to XK 
and \Q would imply T UU (XQ,\K) = \- 2 T uu (Q,K) |27| . 
This scaling could explain this drop of amplitude. How- 
ever the inertial range in our DNS is too small to test 
self-similarity and a large part of the drop is due to the 
presence of viscosity. 

The forced mode has a slightly different behavior. The 
transfer rate from the forced wave number to its nearby 
shells has a considerably larger amplitude. Also, for 
both flows there is some backscattering from the forced 
wave number to shells with smaller wavenumber. This is 
clearer in the helical (ABC) flow. 



Locality of interactions in hydrodynamic turbulence 
have been investigated before in the literature 0, 0, S 
1 rt U^ fl - Although some open issues still remain 
111 Il2 | it has been shown that energy is transfered 
mostly locally. Here, for reasons of comparison we show 
the transfer T UU (Q, K) from hydrodynamical simulations 
using the same external forces and parameters used in 
the MHD simulations. The results are in good agree- 
ment with previous works. 

In figure|21we show the energy transfer for a few modes 
for the TG flow and in figure [21 the energy transfer for 
the ABC flow. In both cases the transfer of energy 
is direct and local: all the curves (with the exception 
of the forced mode Q = 3) are negative for K smaller 
than Q and positive for K larger than Q. As a result, 
all inertial range modes receive energy from modes with 



B. Magneto- Hydrodynamic Turbulence 

We arc now ready to examine results from the en- 
ergy transfer for MHD turbulence. First we examine the 
transfer of kinetic energy from large scales to kinetic en- 
ergy in small scales through the term T UU (Q,K), and 
magnetic energy from large scales to magnetic energy 
in small scales through the term %b(Q,K). These two 
transfer functions bare some significant similarities with 
the hydrodynamic case. 

In figures 0] and [S] we show T uu (top panel) and 7^ 
(bottom panel) for the non-helical TG flow and the heli- 
cal ABC flow. The velocity to velocity transfer has not 
changed drastically (other than a decrease in amplitude) 

UTTaI 



from the pure hydrodynamic case. As in Sec. IIII Al the 
transfer implies a local direct cascade. All the curves are 
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FIG. 4: Top panekThe transfer of energy T UU (Q,K) for the 
Taylor-Green run. The figure shows the rate that kinetic en- 
ergy is transferee! from the modes Q = 3, 10, 20, 30 to kinetic 
energy to all the other modes K. Bottom panel: The trans- 
fer of energy Tbb(Q,K) for the same flow. The figure shows 
the rate that magnetic energy is transfered from the modes 
Q — 3, 10, 20, 30 to magnetic energy to all the other modes 
K. 
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FIG. 5: Same as figure 4 for the ABC run 



negative for K smaller than Q, and positive for K larger 
than Q. Each mode is therefore receiving energy from the 
larger scales (negative transfer) and giving energy to the 
smaller scales (positive transfer). The decrease in am- 
plitude (when compared with the hydrodynamic case) is 
partly because the magnitude of the velocity field is de- 
creased when magnetic field comes to equipartition, and 
partly because now there is a net transfer of energy from 
the velocity field to the magnetic field, making the avail- 
able energy to cascade to small velocity scales smaller. 

The transfer of magnetic energy to magnetic energy 
T bb (Q,K) seems to follow the same behavior as the ve- 
locity field transfer. The results show a direct cascade 
with local transfer of energy from large scales to small 
scales. We note that for the helical case the transfer of 
magnetic energy is larger than the transfer of kinetic en- 
ergy. The likely reason for this behavior is that in the 
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FIG. 6: The transfer of kinetic energy to magnetic energy 
T UU (Q, K) for the Taylor Green run. The figure shows the rate 
that kinetic energy is transfered from the modes Q — 3, 4, 5 
(inset modes Q = 15, 17, 20) to magnetic energy in the modes 
K. 



ABC flow the magnetic energy at large scales and inter- 
mediate scales saturates at higher values than in the TG 
flow, due to the presence of helicity or the dynamo a- 
effect. This process will be discussed in more detail in 
Paper II. 

Next we investigate the transfer of energy from one 
field to the other, by examining the terms T ub and T bu . 
Because of the anti-symmetric property T ub (Q 1 K) = 
—T bu (K,Q), it is sufficient to just study the trans- 
fer of energy from the velocity field to the magnetic 
field. However, we need to remark that unlike the 
1~ UU (Q, K),T bb (Q, K) terms that their dependence on K 
and Q is the same up to a minus sign, the behavior of 
T ub (Q, K) as we vary K is not the same as if we vary Q. 
Therefore the two behaviors need to be studied separately 
(i.e., the transfer of energy from a velocity mode to two 
different magnetic modes is different from the transfer 
of energy from two different velocity modes to a mag- 
netic mode). In figure HH (TG), and [7| (ABC), we show 
the transfer of kinetic energy from the velocity modes 
Q = 3, 4, 5, 15, 17, and 20 to all the examined magnetic 
modes K. 

A few things should be noted. First, in both runs 
(ABC and TG) the modes associated with the large scale 
flow [Q = 3) seem to play a dominant role in the transfer 
of energy from the velocity field to the magnetic field. 
Note also that there is a wider range of magnetic field 
modes into which the forced velocity field modes input 
energy. 

This is more apparent for the helical flow, that seems 
better at stretching and folding the magnetic field. The 
cascade in the modes inside the inertial range is direct in 
both cases but with a small difference. In both cases the 
large scale velocity field is transferring energy to smaller 
scale magnetic field and receiving energy from larger scale 
magnetic field. However, for the Taylor-Green case there 
is very small transfer from one field to the other in the 
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FIG. 7: The transfer of kinetic energy to magnetic energy 
T u b(Q,K) for the ABC run. The figure shows the rate that 
kinetic energy is transferee! from the modes Q — 3, 4, 5 (inset 
modes Q = 15, 17, 20) to magnetic energy in the modes K. 
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FIG. 8: The transfer of kinetic energy to magnetic energy 
T u b(Q, K,) for the Taylor Green run. The figure shows the 
rate that kinetic energy is transfered from the modes Q (x- 
axis) to magnetic energy in the modes Q — 10 (top panel), 
Q — 20 (middle panel), Q = 30 (bottom panel). 



same shell. On the other hand, in the ABC flow the peak 
of the transfer from the magnetic field to the velocity field 
(the negative peaks in figurc[7|) is for the same shell. Note 
also that for the if -shells larger than Q, the transfer for 
all Q follows the same curve. This implies that all the 
small scale velocity modes give energy to the magnetic 
field modes at the same rate. This is clearer when we 
examine the dependence with Q. 

In figures El and El wc show the same transfer function 
T ub (Q, K) for three values of K = 10, 20, 30. The energy 
cascade is also direct (energy going from large scales to 
small scales), however it is clear from these figures that 
the transfer from the velocity field to the magnetic field 
is a highly non-local process. Each magnetic field mode 
Q is receiving energy (positive from all the veloc- 
ity modes with wave number K smaller than Q, with 
the same rate! The only exception is the mechanically 
sustained large scale velocity field that gives even more 








0.0004 
S 0.0000 
-0.0004 

( 

0.0002 " 
Si 0.0000 ■ 
-0.0002 . 







FIG. 9: The transfer of kinetic energy to magnetic energy 
T u t(Q, K, ) for the ABC run. The figure shows the rate that 
kinetic energy is transfered from the modes Q (x-axis) to mag- 
netic energy in the modes Q = 10 (top panel), Q = 20 (middle 
panel), Q = 30 (bottom panel). 



energy (observe the peak at k = 3). In fact, most of the 
energy that is transfered from the velocity field to the 
magnetic field originates from the velocity field modes at 
Q = 3 (around 60% for the TG run and 75% for the ABC 
run.) This energy turns into magnetic energy at several 
wavenumbers K which locally cascades to smaller scales 
through the %{, term. This bigger contribution of the 
large scale flow to 7" u f, (compared with the contribution 
of the turbulent components) is in good agreement with 
the suppression of small scale velocity fluctuations by the 
large scale magnetic field, as observed in [2^. However, 
we need to note that as the scale of the magnetic field be- 
comes smaller there is more energy input from the turbu- 
lent components of the velocity field than from the large 
scale (forced) flow. This just follows from the fact that 
for K large enough, the area below the curve with con- 
stant Tub is larger than the peak at Q — 3. It is possible 
therefore that in the limit of large inertial range the effect 
of the forced velocity scales in the small magnetic scales 
will not be as strong. Finally we note that this mecha- 
nism described above is different in a kinematic dynamo 
regime, as is shown in Paper II. 

In summary, the existence of the long plateau with 
constant T u b(Q,K) at each fixed value of K, and the 
fact that all the magnetic wavenumbers K receive energy 
from the large scale flow at Q = 3 points that interactions 
between the velocity field and the magnetic field are non- 
local in Fourier space. 

This non-local behavior of energy transfer from the ve- 
locity field to the magnetic field seems to be absent from 
the decaying MHD turbulence case studied by J24|. In 
that case although the T u b and 7;, u were more non-local 
than the %b and T uu terms (since energy was transfered 
from the former ones in a wider range of shells than the 
later ones), eventually at large separation of wave num- 
bers the transfer goes to zero. This is very different from 
the plateau behavior wc observe in the forced turbulence 
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FIG. 10: A comparison of the transfers 

T uu (Q,K),T bb ,{Q,K)T ub {Q,K) and T bu (Q,K) for Q = 15 
for the Taylor Green flow. 
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FIG. 11: A comparison of the transfers 

T uu (Q,K),T bb ,(Q,K),T ub (Q,K) and T bu (Q, K) for Q = 15 
for the ABC flow. 

runs. Wc suspect that this difference is due to the fact 
that in the mechanically forced turbulence there is a net 
flux of energy from the velocity field to the magnetic field 
that is responsible for the formation of the plateau which 
does not exist in the decaying turbulence case. 




K, Q=10 

FIG. 12: A comparison of the large K tails of transfers 
T uu (Q,K),T bb (Q,K) and T ub (Q,K) in a log-linear plot for 
Q = 10 for the Taylor Green flow flow. 

u. In the case of the ABC flow the magnitude of the b 
to b transfer seems to be twice the magnitude of the u 
to u transfer. This is due to the fact that the magnetic 
energy in this run is larger than in the TG run at large 
and intermediate scales. 

Figure ^] illustrates the transfer functions 
T uu (Q,K),T bb ,(Q,K) and T ub ,(Q,K) as in figure 
ITUl (TG flow), but we focus here on the large K tail of 
the transfer and we consider Q=10. The fastest drop is 
for the transfer T UU (Q, K) making it the most 'local' one, 
next come the T bb (Q,K) transfer, and finally T ub (Q,K) 
has the slowest drop. The same result was obtained for 
the ABC flow (not shown here). 

Figures ITUl 1111 and 1 121 (as well as a comparison of the 
nonlocal transfers shown in figures [S] and [U] with the local 
transfers in figures 01 and respectively) show that local 
interactions between the same fields are much stronger 
than nonlocal interactions between different fields. How- 
ever, nonlocal interactions spread over several shells, and 
the magnetic field at a given scale K can receive (give) 
energy from (to) several velocity field Q wavenumbers 
(instead of mostly the nearest neighbors as is the case 
for local interactions). Figure ITTjl shows the ratio 



C. A comparison between the transfers 

In the previous section we showed that the transfer 
of energy from the velocity field to velocity field and 
from the magnetic field to magnetic field exhibit a lo- 
cal behavior similar to the transfer in hydrodynamic tur- 
bulence, and the transfer from one field to the other is 
exhibiting a non-local behavior. In order to draw conclu- 
sions we need to compare the magnitude of these trans- 
fers. Figures ITUl and ITU show a comparison of the trans- 
fers T au (Q,K),T bb (Q,K),T ub (Q,K) and T bu (Q,K) with 
Q = 15 for the TG and ABC runs respectively. The 
local transfers u to u and b to b appear to be of larger 
magnitude than the non-local transfers u to b and b to 



— (JT) = T ^(Q, K) / Y, r bb(Q, K). (20) 

Q=l ' Q=l 

This is the ratio of the total energy that the mag- 
netic field at the shell K receives from the velocity field 
through non-local transfer, to the total magnetic energy 
received at the same scale through the local direct cas- 
cade of (magnetic) energy. Although in individual shells 
the local interactions are one order of magnitude larger 
than the non-local transfer, the net amount of energy re- 
ceived at a given scale K by the two processes is compara- 
ble (this ratio is different in a kinematic dynamo regime, 
as will be shown in Paper II). At small scales, the ratio 
seems to settle to a value close to 0.2, indicating that 
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FIG. 13: The ratio NL/L of energy received through the 
non-local transfer T u \, to local 76b, for the ABC and TG sim- 
ulations. The small scales receive 20% of their energy through 
the non-local transfer T u f,. 

20% of the energy received by these scales is through the 
non-local transfer Tub- 



IV. CONCLUSIONS 

In this paper we examined the transfer of energy in 
forced MHD turbulence between the different scales and 
fields involved using the results from numerical simula- 
tions in a turbulent steady state sustained by a mechani- 
cal external force. No qualitative differences in the trans- 
fer of kinetic to kinetic (or magnetic to magnetic) energy 
has been observed, when compared against the transfer 
of energy in a hydrodynamic simulation. These transfers 
were found to be always local and direct. However, all 
kinetic energy modes have been observed to give energy 
to magnetic modes non-locally, in the sense that a small 
scale magnetic field receives the same amount of energy 
from all larger scales of the velocity field in the inertial 
range. Also each magnetic mode was found to receive a 
significant amount of energy from the large scale flow at 
\kp\ ~ 3 (the scale of the forcing), an effect that seems 
to become smaller as we move to smaller scales in the 
inertial range. We note that it is the non-local inter- 
actions that actually sustain the magnetic field against 
Ohmic dissipation. A summary of our results is sketched 
in figure ITU 

We have already noted that a different behavior for 
the non-local transfers 7~ u b and %, u was obtained for the 
mechanically forced turbulence investigated in this work, 
when compared with the decaying turbulence case stud- 
ied by j24|. Compared with incompressible hydrody- 
namic turbulence, involving only one field and one trans- 
fer function, MHD turbulence is richer and more com- 
plex. It involves two interacting fields, several transfer 
functions, and as a result the energy injected at large 
scales can travel to small scales through several channels. 




FIG. 14: A sketch of the energy transfer between different 
scales and different fields. The thickness of the lines is an 
indication of the magnitude of the transfers. The figure il- 
lustrates how energy is transfered to magnetic modes with 
wavenumber k — q in the inertial range. The transfers be- 
tween same fields is always local and direct. Each magnetic 
mode receives energy from all larger in scale velocity modes 
and gives to slightly smaller in scale velocity modes. 
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FIG. 15: A comparison of the Elsasser energy transfer 
T z + Z + (Q, K) with the transfers from the local T uu + %,b and 
the non-local T u b + Tt u contributions. The inset is showing a 
blow up of the small K tail where the non-local interactions 
are more dominant. 

Also the number of quadratic ideal invariants is larger, 
and inverse cascades (not present in three dimensional 
hydrodynamics) can take place. This suggests that in 
MHD flows the particular way the system is set-up (e.g. 
mechanically or magnetically forced, free decaying cases 
without external forces), or even the scale at which the 
energy is injected (compared with the length of the box), 
might have a direct effect in the evolution of the flow and 
lead to different transfers. 

We would also like to comment on the implica- 
tions of our results to the different models of magncto- 
hydrodynamic turbulence. In the pre sent phenomeno- 
logical models of MHD-turbulence Hi Q Hi E E3 , 
locality of the energy transfer is assumed. That is to 
say, these models are derived assuming that scales of dif- 
ferent magnitude do not strongly interact. While this 
assumption seems to be valid for HD-turbulence, this is 
not necessarily true for MHD. As we have shown non- 
local interactions are present in MHD turbulence and 
control the u to b transfers of energy. However, these 
non-local interactions are smaller in amplitude and most 
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of the input of energy to the magnetic field comes from 
the large scale flow and then cascades to smaller scales 
making the assumption of locality justified to some ex- 
tend. However non-local u — to — b need to be considered 
to have a proper description of the energy cascade. To 
illustrate this we show in figure H5I the energy transfer in 
terms of the Elsasscr variables = u ± b. often used 
in turbulence models and we compare it with the contri- 
butions from u to u, b to b, b to u, and u to b. In the 
figure we plot T z + Z + = — / z + Vz + Q<ix and compare 
it with the energy transfer due to the local transfer terms 
%m + %>b an d the non-local transfer terms T u \, + % u . The 
local transfer terms appear to be dominant, except in the 
tails where the transfer of Elsasser variables is dominated 
by the non-local transfers between the magnetic and ki- 
netic energies. This tails, although with small amplitude, 
cannot be completely neglected, as shown by the NL/L 
ratio of figure 13. The non-local tail in the transfer gives 
a net contribution of energy at magnetic small scales of 
roughly 1/5 when compared with the local transfer. 

Finally, we would like to say that our results were based 
on numerical simulations of moderate Reynolds number 



much smaller than what is observed in most physical phe- 
nomena. We already noted that due to the small inertial 
range we cannot test self similarity that would require 
to compare the transfers (i.e. T UU (Q, Ki), T UU {Q,K2) ) 
to wave numbers that are both significantly away from 
each other (K\ -C K2) and away from the forced and 
dissipative scales (Kp <C K\ and K2 <C k v ). 

Finally the transfer of magnetic hclicity and cross- 
hclicity and their effect on the turbulence dynamics is 
also worth studying, but we leave these issues however 
for our future work. 
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